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' Abstract. We establish long-time stability of multi-dimensional noncharacteristic bound- 

OiJ| ary layers of a class of hyperbolic-parabolic systems including the compressible Navier- 

^ ■ Stokes equations with inflow [outflow] boundary conditions, under the assumption of strong 

■^^ ' spectral, or uniform Evans, stability. Evans stabiity has been verified for small-amplitude 

layers by Gues, Metivier, Williams, and Zumbrun. For large-amplitude layers, it may 

be efficiently checked numerically, as done in the one-dimensional case by Costanzino, 

Humpherys, Nguyen, and Zumbrun. 
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2 t. nguyen and k. zumbrun 

1. Introduction 
We consider a boundary layer, or stationary solution, 

(1.1) U = U{xi), lim U{z) = U+, U{0) = Uq 

of a system of conservation laws on the quarter-space 

(1.2) Ut + Y.^'^^U = IZ(^'''(^)^-^)-.' ^ G < = {^1 > 0}, t > 0, 

j jk 

U,F^ G M", B^'' e M"""", with initial data U{x,0) = Uo{x) and Dirichlet type boundary 



conditions specified in (jl.Sp . (|1.6p below. A fundamental question connected to the physical 
motivations from aerodynamics is whether or not such boundary layer solutions are stable 
in the sense of PDE, i.e., whether or not a sufficiently small perturbation of U remains 
close to U, or converges time-asymptotically to U, under the evolution of ()1.2p . That is the 
question we address here. 

1.1. Equations and assumptions. We consider the general hyperbolic-parabolic system 
of conservation laws (II. 2p in conserved variable U, with 

u\ ^ / 



", and {) G M'', where 

jk 

Following |MaZ41 [Z3 l IZ4j. we assume that equations ()1.2p can be written, alternatively, 
after a triangular change of coordinates 

(1.3) W := W{U) = ( j//}^l. 

in the quasilinear, partially symmetric hyperbolic-parabolic form, 

(1.4) A'^Wt + Y^ Am,^ = Yi&'^w,,)^^ + G, 

j jk 

where, defining W+ := W{U+), 

(Al) A^W+),A°,A\^ are symmetric, A° block diagonal, A^ > 9o > 0, 

(A2) for each ^ G M"^ \ {0}, no eigenvector of Y^j S,jA^A^)~^{W+) lies in the kernel of 

Ejk^j^k&HA'^)-Hw+), 



(A3) B^'^ = (^Q ^.,J , E ^%-6 > 0\e, and G = (^^J with giW^, W,) = 0{\W, 
Along with the above structural assumptions, we make the following technical hypotheses: 

(HO) F^,B^'',A'^,A^,&'',W {■),(){■,■) e C^ with s > [{d - l)/2] + 5 in our analysis of 
linearized stability, and s > s{d) := [{d — l)/2] + 7 in our analysis of nonlinear stability. 



STABILITY OF MULTI-DIMENSIONAL BOUNDARY LAYERS 3 

(HI) A\^ is either strictly positive or strictly negative, that is, either A\^ > ^i > 0, or 
A\^ < —6i < 0. (We shall call these cases the inflow case or outflow case, correspondingly.) 

(H2) The eigenvalues of dF^{U^) are distinct and nonzero. 

(H3) The eigenvalues of ^ • dF^^j have constant multiplicity with respect to ^ G M , 

(H4) The set of branch points of the eigenvalues of {A^)~^{iTA^ + J2j^i ^Ci^"')±i t" S K, 
^ G R'^"^ is the (possibly intersecting) union of finitely many smooth curves r = ^g (C), on 
which the branching eigenvalue has constant multiplicity Sq (by definition > 2). 

Condition (HI) corresponds to hyperbolic-parabolic noncharacteristicity, while (H2) is the 
condition for the hyperbolicity at U+ of the associated first-order hyperbolic system obtained 
by dropping second-order terms. The assumptions (A1)-(A3) and (H0)-(H2) are satisfied 
for gas dynamics and MHD with van der Waals equation of state under inflow or outflow 
conditions; see discussions in |MaZ4[ ICHNZl KlMWZSl TGMWZG ]. Condition (H3) holds 
always for gas dynamics, but fails always for MHD in dimension d>2. Condition (H4) is a 
technical requirement of the analysis introduced in [Z2| . It is satisfied always in dimension 
fi = 2 or for rotationally invariant systems in dimensions d >2, for which it serves only to 
define notation; in particular, it holds always for gas dynamics. 

We also assume: 

(B) Dirichlet boundary conditions in PF-coordinates: 

(1.5) iw^,w^^){0,i,t) = h{x,t) := {hiM){5:,t) 
for the inflow case, and 

(1.6) w^\{),x,t) = h{S:,t) 
for the outflow case, with x = {xi,x) G W^. 

This is sufficient for the main physical applications; the situation of more general, Neu- 
mann and mixed-type boundary conditions on the parabolic variable v can be treated as 
discussed in |GMWZ5l[G"MWZ6j . 

Example 1.1. The main example we have in mind consists of laminar solutions {p, u, e){xi,t) 
of the compressible Navier-Stokes equations 

' dtp + div(/9u) = 

dt{pu) + dw{pu^u) + Vp = efiAu + e{p + 77)Vdivn 

dt{pE) + dw{{pE + p)u) = ekAT + e/idiv((u • V)n) 

-|- e{fi + r/)V(u • divn). 



(1.7) 



X G M'^, on a half-space xi > 0, where p denotes density, u G M'^ velocity, e specific internal 
energy, £^ = 6-1- -^ specific total energy, p = p{p, e) pressure, T = T{p, e) temperature, 
p > and 1 7/ 1 < p first and second coefficients of viscosity, k > the coefficient of heat 
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conduction, and e > (typically small) the reciprocal of the Reynolds number, with no-slip 
suction-type boundary conditions on the velocity, 

Uj{0,X2,... ,Xd) = 0, j ^ 1 and ui{0,X2, ■ ■ ■ ,Xd) = V{x) < 0, 

and prescribed temperature, T{0,X2, ■ ■ ■ ,Xd) = Twaii{x)- Under the standard assumptions 
Pp, Te > 0, this can be seen to satisfy all of the hypotheses (A1)-(A3), (H0)-(H4), (B) in 
the outflow case (|1.6p ; indeed these are satisfied also under much weaker van der Waals gas 
assumptions [MaZ4, Z3l lCHNZllGMWZ5t[GMWZ6| . In particular, boundary-layer solutions 
are of noncharacteristic type, scaling as (p,u,e) = {p,u,e){xi/s), with layer thickness ~ e 
as compared to the ~ ^/£ thickness of the characteristic type found for an impermeable 
boundary. 

This corresponds to the situation of an airfoil with microscopic holes through which 
gas is pumped from the surrounding flow, the microscopic suction imposing a flxed normal 
velocity while the macroscopic surface imposes standard temperature conditions as in flow 
past a (nonporous) plate. This conflguration was suggested by Prandtl and tested experi- 
mentally by G.I. Taylor as a means to reduce drag by stabilizing laminar flow; see [Si IBra] . 
It was implemented in the NASA F-16XL experimental aircraft program in the 1990's with 
reported 25% reduction in drag at supersonic speeds [Braj l I Possible mechanisms for this 
reduction are smaller thickness ~ e << ^/e of noncharacteristic boundary layers as com- 
pared to characteristic type, and greater stability, delaying the transition from laminar to 
turbulent flow. In particular, stability properties appear to be quite important for the 
understanding of this phenomenon. For further discussion, including the related issues 
of matched asymptotic expansion, multi-dimensional effects, and more general boundary 
configurations, see |GMWZ5] . 

Example 1.2. Alternatively, we may consider the compressible Navier-Stokes equations 
()1.7p with blowing-type boundary conditions 



Uj{0,X2,... ,Xd) = 0, j ^ 1 and ui{0,X2, ■ ■ ■ ,Xd) = V{x) > 0, 
and prescribed temperature and pressure 

T{0,X2,...,Xd) =Twau{x), p(0,X2,... , Xrf) = Puiall{x) 

(equivalently, prescribed temperature and density). Under the standard assumptions pp, 
Te > on the equation of state (alternatively, van der Waals gas assumptions), this can be 
seen to satisfy hypotheses (A1)-(A3), (H0)-(II4), (B) in the inflow case (II. Sh . 



Lemma 1.3 i |MaZ3[ [Z3l fCJMWZSl INZ] ^ . Given (A1)-(A3) and (H0)-(H2), a standing wave 
solution (II. ip of (II. 2p . (B) satisfies 



(1.8) {d/dxif{U-U- 



<C7e-^^\ 0<yfc<s + l, 



as xi — > -|-oo, s as in (HO). Moreover, a solution, if it exists, is in the inflow or strictly 
parabolic case unique; in the outflow case it is locally unique. 

Proof. See Lemma 1.3, |NZ| . D 



See also NASA site http://www.dfrc.nasa.gov/Gallery/photo/F-16XL2/index.htnil 
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1.2. The Evans condition and strong spectral stability. The hnearized equations of 
(fr2]l . (B) about U are 

(1-9) Ut = LU := J^iB'^U,,)., " T^i^'U),, 

j,k j 

with initial data U{0) = Uq and boundary conditions in (linearized) VF-coordinates of 

W^(0,5,t) := {w^,w^^)'^{<d,S:,t) = h 

for the inflow case, and 

w^\Q,5:,t) = h 

for the outflow case, with x = (xi,x) G M'^, where W := {dW / dU){tJ)U . 

A necessary condition for linearized stability is weak spectral stability, defined as nonexis- 
tence of unstable spectra 5RA > of the linearized operator L about the wave. As described 
in Section [2.1.11 this is equivalent to nonvanishing for all ^ G M , 3f?A > of the Evans 
function 

(defined in (j2.8p ). a Wronskian associated with the Fourier-transformed eigenvalue ODE. 

Definition 1.4. We define strong spectral stability as uniform Evans stability: 

(D) \DL{i,X)\>eiC)>0 

for {i, A) on bounded subsets C C {| G M"'~\ KA > 0} \ {0}. 

For the class of equations we consider, this is equivalent to the uniform Evans condition of 
[GMWZSl IGA/IWZ6] , which includes an additional high-frequency condition that for these 
equations is always satisfied (see Proposition 3.8, |GMWZ5] ). A fundamental result proved 
in |GMWZ5] is that small-amplitude noncharacteristic boundary-layers are always strongly 
spectrally stableo 

Proposition 1.5 r |GMWZ5j ). Assuming (A1)-(A3), (H0)-(H3), (B) for some fixed end- 
state (or compact set of endstates) C/+, boundary layers with amplitude 

\W ~ ^+IIl°°[0,+oo] 

sufficiently small satisfy the strong spectral stability condition (D). 



As demonstrated in [SZ| . stability of large-amplitude boundary layers may fail for the 
class of equations considered here, even in a single space dimension, so there is no such 
general theorem in the large-amplitude case. Stability of large-amplitude boundary-layers 
may be checked efficiently by numerical Evans computations as in [BDG[ IBrH IBr21 IBrZ[ 
E^ EHRZI ELZI [CHNZl IHLyZlj rHLyZ2| . 



The result of [GMWZS] applies also to more general types of boundary conditions and in some situations 
to systems with variable multiplicity characteristics, including, in some parameter ranges, MHD. 
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1.3. Main results. Our main results are as follows. 

Theorem 1.6 (Linearized stability). Assuming (A1)-(A3), (H0)-(H4), (B), and strong 
spectral stability (D), we obtain asymptotic L^ n //K"^^)/2|+^ -^ LP stability of (jl.Op in 
dimension d >2, and any 2 < p < oo, with rate of decay 

\U{t)\L2 < C(l + t)-^(|C/o|iinH3 +i?o), 

\U{t)\L. < Cil + t)--2(^-yP^+y'P{\UoynHii,~i)f2n^. + Eo), 

provided that the initial perturbations Uq are in L} n H^ for p = 2, or in L^ n //I('^^i)/21+5 
for p > 2, and boundary perturbations h satisfy 

\hit)\L2<Eoii+t)-^''+'y\ 

X 

(1.11) \h{t)\L^<Eo{l + t)-''/^ 

|2?h(t)LinH[(d-i)/2i+5 < Eoil + t)-''/^-', 

X X 

where T>fi{t) := |/it| + |/ii| + \hxx\, Eq is some positive constant, and e > is arbitrary small 
for the case d = 2 and e = for d> 3. 

Theorem 1.7 (Nonlinear stability). Assuming (Al)-(A3), (H0)-(H4), (B), and strong 
spectral stability (D), we obtain asymptotic L} n H^ — > L^ n H"^ stability of U as a solution 
of (II. 2p in dimension d>2, for s > s{d) as defined in (HO), and any 2 < p < oo, with rate 
of decay 

\m - U\l. < Cil + t)-i(l-l/rt+l/2p(|[;^|^,^^^ + Eo) 
\U{t) - U\hs < C(l + tr^{\Uo\LinHs + Eo), 

provided that the initial perturbations Uo '■= Uo ~ U are sufficiently small in L} Pi H^ and 
boundary perturbations h{t) := h{t) — W{Uo) satisfy (jl.lip and 



d-l 
4 



(1.13) Bh{t)<Eo{l+t) 

with sufficiently small Eo, where the boundary measure B^ is defined as 

[(«+l)/2] 

(1-14) BH{t):=\h\Hs^s:)+ Yl \91^\lH^) 

for the outflow case, and similarly 

[(s+l)/2] 

(1.15) Bhit):=\h\Hs(i)+ Yl l^i^2|L2(x)+Xll^t^il^'« 

i=0 i=0 

for the inflow case. 

Combining Theorem 11.71 and Proposition 11.51 we obtain the following small-amplitude 
stability result, applying in particular to the motivating situation of Example ll.il 
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Corollary 1.8. Assuming (A1)-(A3), (H0)-(H4), (B) for some fixed endstate (or compact 
set of endstates) C/+, boundary layers with amplitude 

W ~ ^+IIl°°[0,+oo] 

sufficiently small are linearly and nonlinearly stable in the sense of Theorems \1.6\ and \1.7\ 

Remark 1.9. The obtained rate of decay in L^ may be recognized as that of a {d — 1)- 
dimensional heat kernel, and the obtained rate of decay in L°° as that of a d-dimensional 
heat kernel. We believe that the sharp rate of decay in L^ is rather that of a d-dimensional 
heat kernel and the sharp rate of decay in L°° dependent on the characteristic structure of 
the associated inviscid equations, as in the constant-coefficient case |HoZH lHoZ2j . 

Remark 1.10. In one dimension, strong spectral stability is necessary for linearized as- 
ymptotic stability; see Theorem 1.6, ^NZj . However, in multi-dimensions, it appears likely 
that, as in the shock case |Z3j , there are intermediate possibilities between strong and weak 
spectral stability for which linearized stability might hold with degraded rates of decay. In 
any case, the gap between the necessary weak spectral and the sufficient strong spectral 
stability conditions concerns only pure imaginary spectra JRA = on the boundary between 
strictly stable and unstable half-planes, so this should not interfere with investigation of 
physical stability regions. 

1.4. Discussion and open problems. Asymptotic stability, without rates of decay, has 
been shown for small amplitude noncharacteristic "normal" boundary layers of the isen- 
tropic compressible Navier-Stokes equations with outflow boundary conditions and vanish- 
ing transverse velocity in |KK| . using energy estimates. Corollary 11.81 recovers this existing 
result and extends it to the general arbitrary transverse velocity, outflow or inflow, and 
isentropic or nonisentropic (full compressible Navier-Stokes) case, in addition giving as- 
ymptotic rates of decay. Moreover, we treat perturbations of boundary as well as initial 
data, as previous time- asymptotic investigations (with the exception of direct predecessors 
[YZtlNZ] ) do not. As discussed in Appendix [Aj the type of boundary layer relevant to the 
drag-reduction strategy discussed in Examples I1.1H1.2I is a noncharacteristic "transverse" 
type with constant normal velocity, complementary to the normal type considered in JKKJ. 

The large-amplitude asymptotic stability result of Theorem 11.71 extends to multi dimen- 
sions corresponding one-dimensional results of |YZ[INZJ . reducing the problem of stability to 
verification of a numerically checkable Evans condition. See also the related, but technically 
rather different, work on the small viscosity limit in [ MZl IGMWZSj IGMWZ6J . By a combi- 
nation of numerical Evans function computations and asymptotic ODE estimates, spectral 
stability has been checked for arbitrary amplitude noncharacteristic boundary layers of the 
one-dimensional isentropic compressible Navier-Stokes equations in |CHNZ| . Extensions to 
the nonisentropic and multi-dimensional case should be possible by the methods used in 
|HLyZl| and |HLyZ2| respectively to treat the related shock stability problem. 

This (investigation of large-amplitude spectral stability) would be a very interesting di- 
rection for further investigation. In particular, note that it is large-amplitude stability that 
is relevant to drag-reduction at flight speeds, since the transverse relative velocity (i.e., 
velocity parallel to the airfoil) is zero at the wing surface and flight speed outside a thin 
boundary layer, so that variation across the boundary layer is substantial. We discuss this 
problem further in Appendix [X] for the model isentropic case. 
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Our method of analysis follows the basic approach introduced in |Z21 IZ31 IZ4j for the 
study of multi-dimensional shock stability and we are able to make use of much of that 
analysis without modification. However, there are some new difficulties to be overcome in 
the boundary-layer case. 

The main new difficulty is that the boundary-layer case is analogous to the undercom- 
pressive shock case rather than the more favorable Lax shock case emphasized in |Z3] . in 
that Gy-^ 9^ t~^/'^G as in the Lax shock case but rather Gy^ ~ (e~^l^il +t~^/'^)G^ ^ > 0, as in 
the undercompressive case. This is a significant difficulty; indeed, for this reason, the un- 
dercompressive shock analysis was carried out in [Z3] only in nonphysical dimensions d > 4. 
On the other hand, there is no translational invariance in the boundary layer problem, so 
no zero-eigenvalue and no pole of the resolvent kernel at the origin for the one-dimensional 
operator, and in this sense G is somewhat better in the boundary layer than in the shock 
case. 

Thus, the difficulty of the present problem is roughly intermediate to that of the Lax 
and undercompressive shock cases. Though the undercompressive shock case is still open in 
multi-dimensions for c? < 3, the slight advantage afforded by lack of pole terms allows us to 
close the argument in the boundary-layer case. Specifically, thanks to the absence of pole 
terms, we are able to get a slightly improved rate of decay in L°°(a;i) norms, though our 
LP'{xi) estimates remain the same as in the shock case. By keeping track of these improved 
sup norm bounds throughout the proof, we are able to close the argument without using 
detailed pointwise bounds as in the one-dimensional analyses of [HZtlRZ) . 

Other difficulties include the appearance of boundary terms in integrations by parts, 
which makes the auxiliary energy estimates by which we control high-frequency effects 
considerably more difficult in the boundary-layer than in the shock-layer case, and the 
treatment of boundary perturbations. In terms of the homogeneous Green function G, 
boundary perturbations lead by a standard duality argument to contributions consisting of 
integrals on the boundary of perturbations against various derivatives of G, and these are 
a bit too singular as time goes to zero to be absolutely integrable. Following the strategy 
introduced in [YZ[ INZJ . we instead use duality to convert these to less singular integrals 
over the whole space, that are absolutely integrable in time. However, we make a key 
improvement here over the treatment in |YZ1 INZj . integrating against an exponentially 
decaying test function to obtain terms of exactly the same form already treated for the 
homogeneous problem. This is necessary for us in the multi-dimensional case, for which 
we have insufficient information about individual parts of the solution operator to estimate 
them separately as in |YZ1 INZ] , but makes things much more transparent also in the one- 
dimensional case. 

Among physical systems, our hypotheses appear to apply to and essentially only to the 
case of compressible Navier-Stokes equations with inflow or outflow boundary conditions. 
However, the method of analysis should apply, with suitable modifications, to more general 
situations such as MHD; see for example the recent results on the related small-viscosity 
problem in [GMWZSJ IGMWZ6] . The extension to MHD is a very interesting open problem. 

Finally, as pointed out in Remark 11.101 the strong spectral stability condition does not 
appear to be necessary for asymptotic stability. It would be interesting to develop a refined 
stability condition similarly as was done in |SZ1 \7j2\ IZ3t IZ4| for the shock case. 
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2. Resolvent kernel: construction and low-frequency bounds 

In this section, we briefly recal the construction of resolvent kernel and then establish 
the pointwise low-frequency bounds on G? ^, by appropriately modifying the proof in Z3] 
in the boundary layer context |YZ1 INZj . 

2.1. Construction. We construct a representation for the family of elliptic Green distri- 
butions G;^(xi,yi), 

(2.1) G^,{;yi):={L^-X)-'6y,{.), 

associated with the ordinary differential operators (L; — A), i.e. the resolvent kernel of 
the Fourier transform Lx of the linearized operator L of ()1.9p . To do so, we study the 
homogeneous eigenvalue equation (L^ — A)f7 = 0, or 

LoU 

(2.2) m m 

+ iJ2iB"'CkUy - ^ B^%CkU -XU = 0, 
with boundary conditions (translated from those in 1^-coordinates) 

(-) (""-T"'"'''4><°>^(o 

where * = for the inflow case and is arbitrary for the outflow case. 
Define 

n 

A^ := n ^n^) 

i=i 
where A^(^) denote the open sets bounded on the left by the algebraic curves A^(.^i,^) 
determined by the eigenvalues of the symbols —^^B^ — i^A^ of the limiting constant- 
coefficient operators 

Lx w := B^w" — A^ui' 

as ^1 is varied along the real axis, with ^ held fixed. The curves A^(-,^) comprise the 
essential spectrum of operators ^£_,_- Let A denote the set of (^, A) such that A G A^. 
For (^,A) G A^, introduce locally analytically chosen (in ^, A) matrices 

(2.4) <!>+ = (</,+ ,...,</>+), $o = (</,o+i,...,</>o+,), 
and 

(2.5) $ = ($+,$°), 

whose columns span the subspaces of solutions of ()2.2p that, respectively, decay at x = -|-oo 
and satisfy the prescribed boundary conditions at x = 0, and locally analytically chosen 
matrices 

(2.6) ^i,^ = (^<l,...,^l), $+ = (V.+ ,...,V^++,) 
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and 

(2.7) ^ = (^°,^+). 

whose columns span complementary subspaces. The existence of such matrices is guaranteed 



by the general Evans function framework of |AGJ1 IGZl IMaZ3| ; see in particular |Z3^ INZI 
That dimensions sum to n + r follows by a general result of [GMWZ5] ; see also 



2.1.1. The Evans function. Following [AGJllGZl ISZj . we define on A the Evans function 

(2.8) i?L(|,A):=det($o,$+)|,.=o. 

Evidently, eigenfunctions decaying at +oo and satisfying the prescribed boundary conditions 
at xi = occur precisely when the subspaces span*!?*^ and span*!'"'" intersect, i.e., at zeros 
of the Evans function 

^l(I,A) = 0. 

The Evans function as constructed here is locally analytic in (^, A), which is all that we 
need for our analysis; we prescribe different versions of the Evans function as needed on 
different neighborhoods of A. Note that A includes all of {| G IR'^~^ SRA > 0} \ {0}, so that 
Definition 11.41 is well-defined and equivalent to simple nonvanishing, away from the origin 
(^, A) = (0,0). To make sense of this definition near the origin, we must insist that the 
matrices $-^ in ()2.8I) remain uniformly bounded, a condition that can always be achieved by 
limiting the neighborhood of definition. 

For the class of equations we consider, the Evans function may in fact be extended 
continuously along rays through the origin [R2 | IMZj [GMWZSj [GMWZ6] . 

2.1.2. Basic representation formulae. Define the solution operator from yi to xi of ODE 
(L| - X)U = 0, denoted by J^v^^^^ as 

^J/i--i = cI>(xi,A)$-i(yi,A) 

and the projections 11^^,11^^ on the stable manifolds at 0,+oo as 

n+=(<5+(yi) 0)<l>-Hyi), < = (0 <^^{yi))<l>~Hyi). 

We define also the dual subspaces of solutions of {L*- — \*)W = 0. We denote growing 
solutions 

(2.9) ^o = (<^;,---,0O), ^+ = {4>t+,r--,4>t+r), 
$ := ($'^, ^"'") and decaying solutions 

(2.10) §o^(v;?,...,^+), §+ = (v5+,i,---,v5++.), 

and 4/ := {^^, ^^), satisfying the relations 

(2.11) (* ^);jH^ ^)o,+ = ^' 

where 

(2.12) S^ = 
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With these preparations, the construction of the Resolvent kernel goes exactly as in the 
construction performed in [ZH^ IMaZ3l IZ3| on the whole line and [YZ[ INZ| on the half line, 
yielding the following basic representation formulae; for a proof, see |MaZ3tlNZ| . 

Proposition 2.1. We have the following representation 

(^,^. ^ t /(/„,, 0).F^l--l^+(5«r^yl)Un,0)*^ for x^>y^, 

^Z.L6) U^^,[x^,y^) j_(^^^o)_^,,^x,^o^(5C)-l(y,)(J„,0)*^ for x,<y,. 

Proposition 2.2. The resolvent kernel may alternatively he expressed as 

f(/„,0)cl>+(xi;A)M+(A)^o*(yi;A)(4,Or^ xi > yi, 
«'^^^''^'^ \-(/„,0)cl>0(xi;A)MO(A)^+*(yi;A)(/„,Or xi < yi, 
where 

(2.14) M(A) := diag(M+(A),M°(A)) = ^-^{z]X){S^-^{z)i>-^*{z]X). 

2.1.3. Scattering decomposition. From Propositions 12.11 and 12.21 we obtain the following 
scattering decomposition, generalizing the Fourier transform representation in the constant- 
coefficient case, from which we will obtain pointwise bounds in the low-frequency regime. 

Corollary 2.3. On A^ n p{L^), 

(2.15) G^^^{xuyi) = Y.dl,(t>t{xi;X)i'tiyi;Xr + Y.<Pti^i;>^)4'tiyi-^>^T 

j,k k 

for < yi < xi, and 

(2.16) G^-,(xi,yi) = J^dO,0+(xi;A)^^+(yi;A)*-5^Vfc+(^i;A)^^+(yi;A)* 

j,k k 

for < xi < yi, where 

(2.17) d°^+(A) = (/,0)($+ $0)"'^+. 
Proof. For < xi < yi, we obtain the preliminary representation 

Gi,{xuyi) = j;4(A)<A+(xi;A)v;+(yi;A)* + j;eO,v;(xi;A)v;+(yi;A)* 

j,k jk 

from which, together with duality (j2.1ip . representation (j2.13p . and the fact that Ho = 
/ — n+ , we have 



(2.18) =-($+ ^+) ^ 



1 = - (^1.+ §+) ^no^+ 



/-($+ 0)($+ $0)"^ 



Similarly, for < yi < xi, we obtain the preliminary representation 

G|,(xi,yi) = ^d+(A),/.+ (xi;A)V;+(yi;A)* + 5]e+4(xi;A),^+(yi;A)* 

j,k jk 
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from which, together with duahty (j2.1ip and representation (|2.13p . we have 



(2.19) 



$+Mn+ ("f + $+) 

($+)-! (<5+ o)($+ $o)"^(^+ $+) 

(/ 0)($+ $0)"'(*+ «5+) 



(V o)(^^ ^°)"'^" 



\h 1 



D 



Remark 2.4. In the constant-coefficient case, with a choice of common bases ^'''^ = $+''^ 
at 0,-|-oo, the above representation reduces to the simple formula 

(2 20) G- fxi m) - /^f=A.+i'^.^(^i;^)'^r(yi;^) ^1 > yi' 

^' [-Ej=iV'j (2;i;A)V'7 (yi;A) xi<yi. 

2.2. Pointwise loAv-frequency bounds. We obtain pointwise low-frequency bounds on 
the resolvent kernel G^ ^{xi,yi) by appealing to the detailed analysis of [Z2 1 IZ3 1 [GMWZl] 
in the viscous shock case. Restrict attention to the surface 

(2.21) r«":={A : SReA = -0i(|eP + l^mAp)}, 
for ^1 > sufficiently small. 

Proposition 2.5 (|Z3]). Under the hypotheses of Theorem \1.7\ for A G F^ and p := |(^, A)|, 

9i > 0, and 9 > sufficiently small, there hold: 

(2.22) \G-^^^{xi,yi)\<C^2e-'''"\''-'-y'\. 
and 

(2.23) \dlG^,^{xi,yi)\ < C^2{p^ + l3e''y^)e~'P'\--^-y^\ 
where 

(2.24) ^2 := 1 + j; [p-^|9mA - r/+(0| + p] ^'''''\ 

j 

and Sj,r]'^{^) are as defined in (H4). 

Proof. This follows by a simplified version of the analysis of [Z3J , Section 5 in the viscous 
shock case, replacing $~, ^~ with $'^, ^'^, omitting the refined derivative bounds of Lem- 
mas 5.23 and 5.27 describing special properties of the Lax and overcompressive shock case 
(not relevant here), and setting £ = 0, or 7 = 1 in definition (5.128). Here, i is the mul- 
tiplicity to which the Evans function vanishes at the origin, (^, A) = (0,0), evidently zero 
under assumption (D). The key modes $"*", ^^ at plus spatial infinity are the same for the 
boundary-layer as for the shock case. 

This leads to the pointwise bounds (5.37)^(5.38) given in Proposition 5.10 of [Z3| in case 
a = 1, 71 = 1 corresponding to the uniformly stable undercompressive shock case, but 
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without the first 0{p^^), or "pole", terms appearing on the righthand side, which derive 
from cases 7 '^ p^^ not arising here. But, these are exactly the claimed bounds (j2.22p - 

We omit the (substantial) details of this computation, referring the reader to |Z3j . How- 
ever, the basic idea is, starting with the scattering decomposition of Corollary [2?L3l to note, 
first, that the normal modes $■', ^■', $-^, ^^ can be approximated up to an exponentially 
trivial coordinate change by solutions of the constant-coefficient limiting system at x — > -|-oo 
(the conjugation lemma of |MZj ) and, second, that the coefficients Mj^, djk may be well- 
estimated through formulae (j2.14p and (j2.17p using Kramer's rule and the assumed lower 
bound on th Evans function \D\ appearing in the denominator. This is relatively straight- 
forward away from the branch points 9 A = r]j{^) or "glancing set" of hyperbolic theory; 
the treatment near these points involves some delicate matrix perturbation theory applied 
to the limiting constant-coefficient system at j; — > +00 followed by careful bookkeeping in 
the application of Kramer's rule. D 

3. Linearized estimates 
We next establish estimates on the linearized inhomogeneous problem 

(3.1) Ut-LU = f 

with initial data U{0) = Uq and Dirichlet boundary conditions as usual in T^-coordinates: 

(3.2) W{0, X, t) := {w^,w^^f{0, x,t) = h 
for the inflow case, and 

(3.3) w^^{0,x,t) = h 
for the outflow case, with x = (xi,x) € M*^. 

3.1. Resolvent bounds. Our first step is to estimate solutions of the resolvent equation 
with homogeneous boundary data h = 0. 

Proposition 3.1 (High-frequency bounds). Given (Al)-(A2), (H0)-(H2), and homoge- 
neous boundary conditions (B), for some R, C sufficiently large and 6 > sufficiently small, 

(3-4) I(^|-A)-V1hi(.o^^I/IhM..)' 

and 

for all |(C, A)| > R and 3f?eA > —9, where f is the Fourier transform, of f in variable x and 

\f\mi,,)--=\i^ + \d.A + \i\)f\LH.,)- 

Proof. First observe that a Laplace-Fourier transformed version with respect to variables 
(A, x) of the nonlinear energy estimate in Section llTTl with s = 1, carried out on the linearized 
equations written in ly-coordinates, yields 

(3.6) (5ReA + 0i)|(l + III + \d,,\W\^ < C(\W\^ + (1 + |I|')|H^||/| + \d,,W\\d,J\ 
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for some C big and 6*1 > sufficiently small, where |.| denotes l-lLa^^,^). Applying Young's 
inequality, we obtain 

(3.7) (3f?eA + 0i)|(l + III + \d^,\)W\^ < C\W\^ + C\{1 + ||| + |9.J)/?. 

On the other hand, taking the imaginary part of the L^ inner product of U against 
\U = / + LU, we have also the standard estimate 

(3.8) \QmX\\U\l,<C\U\j,,+C\f\l„ 
and thus, taking the Fourier transform in x, we obtain 

(3.9) \QmX\\W\^ < C\f\^ + C|(l + III + \d,,\)W\''. 
Therefore, taking 6 = 9i/2, we obtain from (j3.7p and (|3.9p 

(3.10) 1(1 + |A|i/2 + III + \d,,\)W\^ < C\W\' + C\{1 + III + |a,J)/T, 

for any 3f?eA > —9. Now take R sufficiently large such that \W\'^ on the right hand side 
of the above can be absorbed into the left hand side, and thus, for all |(C,A)| > R and 
KeA > -9, 

(3.11) 1(1 + |A|^/2 + III ^ |a,j)H^|2 < c|(i + III + \d,M', 

for some large C > 0, which gives the result. D 

We next have the following: 

Proposition 3.2 (Mid-frequency bounds). Given (A1)-(A2), (H0)-(H2), and strong spec- 
tral stability (D), 

(3.12) |(L^--A)"^|j^i(^^) < C7, for R-^ <\{t\)\<R and^e\>-9, 

for any R and C = C{R) sufficiently large and 9 = 9{R) > sufficiently small, where 
\f\fji(^ ) is defined as in Proposition \&.1[ 

Proof. Immediate, by compactness of the set of frequencies under consideration together 
with the fact that the resolvent (A — L;)^^ is analytic with respect to H^ in (^,A); see 



Proposition 4.8, |Z4]. D 

We next obtain the following resolvent bound for low-frequency regions as a direct con- 
sequence of pointwise bounds on the resolvent kernel, obtained in Proposition 12. 5i 

Proposition 3.3 (Low-frequency bounds). Under the hypotheses of Theorem \1.1\ for A G 
r^ and p := |(.^, A)|, 9i sufficiently small, there holds the resolvent bound 

(3.13) |(L^-- A)-iaf jliP(,,) < C72p-'/^[/|/|Li(x,) +/^l/U-(-i)" 

for all 2 < p < cxD, /? = 0, 1, where 72 is as defined in (j2.24p . 
Proof. Using the convolution inequality \g * HIlp < \g\Lp\h\Li and noticing that 
KG,,{x,,y^)\ < C^2{p^ + Pe-'y^)e-'P'\---y-\, 
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we obtain 



(3.14) 



as claimed. 



|(L-- A)-i9fjlip(x,) = / d^^G,Jxi,y,)f{yuOdy 



LP{xi) 



< 



C72(/ + /3e- 



'Syi\„''0p'^\xi~yi\ 



\fiyi,i)\dy, 



LP 



<C^2P-"'' Pn/|Li(xO+/5|/lL-(x.O 



D 



Remark 3.4. The above L^ bounds may alternatively be obtained directly by the argu- 
ment of Section 12, [GMWZT] . using quite different Kreiss symmetrizer techniques, again 
omitting pole terms arising from vanishing of the Evans function at the origin, and also the 
auxiliary problem construction of Section 12.6 used to obtain sharpened bounds in the Lax 
or overcompressive shock case (not relevant here). 

3.2. Estimates on homogeneous solution operators. Define low- and high-frequency 
parts of the linearized solution operator S{t) of the linearized problem with homogeneous 
boundary and forcing data, f,h = 0, as 

1 



(3.15) 



Slit) :-- 



{2my J\i\<rJr^n{\x\<r} 
and 

(3.16) S2{t) := e 

Then we obtain the following: 



^Lt 



e^'+"^-^{L^-Xy'dXd^ 



Slit). 



Proposition 3.5 (Low-frequency estimate). Under the hypotheses of Theorem \1.7\ for 
P= (/3i,/3') with Pi = 0,1, 



|5i(t)5f/|i. <c(i + t)-('^-i)/^-i'5|/Vla + c/3i(i + t)-('^-i)/^|/Ui,o 



(3.17) 



-(d+l)/4-|/3|/2| 



-(d+l)/4| 



|5i(t)af /|^.,oo <c(i + t)-^''+^^/^-i^i/1/|ii + cpiii + t)-^''+'nf\L^^- , 

X.XJ^ x,x^ 



|5i(t)af/|L» <c{i + t] 



-d/2-\!3\/2 



i/iM+c/3i(i+r'/Vh 



where \ ■ \ip,'i denotes the norm in L'^ [x; L'^ [xi)) . 

x,x-^ 

Proof. The proof will follow closely the treatment of the shock case in fZ3]. Let 'u(xi,^, A) 
denote the solution of (L; — \)u = /, where /(xi,^) denotes Fourier transform of /, and 



u{x,t):=Si{t)f 



e^t+^^-(L,-X)-'f{xi,i)dXdl 



C^T^i) J\i\<rJrin{\x\<r} 
Recalling the resolvent estimates in Proposition 13.31 we have 

H^ui^)\LPix,) < C72p-'/^|/|Li(x,) < C72p-'/^|/|Li(x) 

where 72 is as defined in (j2.24p . 
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Therefore, using Parseval's identity, Fubini's theorem, and the triangle inequahty, we 
may estimate 



1 



Hl^xux)(^) - 'J2^ 



XI J £, 



r«n{|A|<r} 



e^%{xi,i,X)d\ 



2 _ 



d^dx\ 



1 



< 



(27r)2'i 

1 
(2^ 



r«n{|A|<r} 



e^^u{xi,i,\)d\ 



L2{xi) 



di 



r«n{|A|<r} 



e^'^'\u{xuiA)\L^ix,)dX 



2 , 
di 



<c|/lii(.) 



e'''^'-f2p-^d\ 



r«n{|A|<r} 



di. 



Specifically, parametrizing F^ by 

A(e,A;) =ifc-ei(A;2 + |eT), A; G M, 
and observing that by ()2.24p . 

|A:-r,-(f)Ki/«.-i 



(3.18) 



i 



|fc-r,(|)|v-i 



where e := \ — - (0 < e < 1 chosen arbitrarily if there are no singularities), we estimate 



r«n{|A|<r} 



e'^'="*72/5-'dA 



2 _ 



di < 



^(>^'+\i\')t^,p-'dk 



2 , 
d^ 



< / g-2ei|5pt|||-2e 



g-^l/c t\J,\e-l^f^ 



d^ 



+ E 



~20i|C|2t|t|-2e 



ler 



g-yife f|^_^^.^_^^|.-i^^ 



2 _ 






e-'^i'^ '\k\'-'dk 



di 



Likewise, we have 



|m|^2,oo {t) 



(2^)2rf 



.At 



r«n{|A|<r} 



e^'n(a;i,^,A)(iA 



L°°{xi) 



dC 



< 



(27r)2rf 



.SReAti -, 



< c^l/li^w 



r«n{|A|<r} 

( 

r«n{|A|<r} 



e |tt(2;i,,^,A)|ioo(^^)(iA 
2 _ 



d^ 



e^^^S2dA 



dC 
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where 



rfn{|A|<r} 



e'^^%rfA'de< /e-2''^l«l'* 



-OikH 



dk 



di 



^?i- 



-2ei|C'pt|||2-2e 



\k-Tj{i)\'-^dk di 



< Ct-^d+l)/2 ^ ^ / g-2ei||pt|||2~2e 



*|A:|^-i(iA; 



de 



Similarly, we estimate 



1 



i^"^^"^.(*)^(2;o^y,- 



„At 



r«n{|A|<r} 



m(xi,^, A)(iA 



L°°(xi) 



di 



< 



Me\t 



{2ttY JiJrin{\x\<r} 



u{xi,i,X)\L^^^^)dXdi 



„KeA< 



< C\f\m,) / f , e'«^^*72C?AdC 

i^Jr«n{|Al<r} 



where as above we have 
J? ir«n{|A|<r} 



e^'^^St^Ad|< /e-^il«'l'* / e'^^^^'dkdi 






Ji JR 



'dkdi 



The xi-derivative bounds follow similarly by using the resolvent bounds in Proposition 

Hwith /?! = 1. The x-derivative bounds are straightforward by the fact that d^f = {iC)^f- 
Finally, each of the above integrals is bounded by C|/|/^i(^.) as the product of \f\Li(^x) 
times the integral quantities 72/0~^, 72 over a bounded domain, hence we may replace t by 
(1 + t) in the above estimates. D 

Next, we obtain estimates on the high-frequency part S2{t) of the linearized solution 
operator. Recall that S2{t) = S{t) — Si{t), where 



S{t) 



{27riy J^, 



S-^e^i'di 



and 



Slit) 



{2'KiY J\S,\<rJnr\{\\\<r} 



e^i+'i-^{L,-X)-^dXdC. 



18 



T. NGUYEN AND K. ZUMBRUN 



Then according to [Z4| Corollary 4.11], we can write 



(3.19) 



S2it)f 



1 /■— Pl+lOO 

P.V, 



X e 



X|?-|2 



1— lOO 



5|2+|SmA|2>6(i+6»2 



Proposition 3.6 (High-frequency estimate). Given (A1)-(A2), (H0)-(H2), (D), and ho- 
mogeneous boundary conditions (B), for < |a| < s — 3, s as in (HO), 

, , |52(i)/U2<Ce-^^*|/b3, 

(3.20) fl , 

|9°52(i)/U2<Ce-^^*|/I^H+3. 

Proof. The proof starts with the following resolvent identity, using analyticity on the resol- 
vent set p{L/) of the resolvent (A — L;)^^, for all / G V{Lc), 

(3.21) (A - L^-)-V = \-\\ - L^r'L^f + A-V. 
Using this identity and (13.19p . we estimate 



S2{t)f 



1 /■— Pl+lOO 

rP.V. 



_i ^|^P+|QmA|2>6(i+e»2 



(3.22) 



(27ri)<= .y-yi-»oo 

X e*«-*+^*A-i(A - L^)-^L^f\xi,^)dCdX 

1 /• — 6*1+400 

+ tt^-t-tP.V. 



(27ri) 



d 



/I— lOO 



_i ^|C|2 + 13'mA|2>6li+6l2 



X e 



j^-x+At \ — 1 



X-'f{xi,Od^dX 



— : S*! -|- 5*2, 
where, by Plancherel's identity and Propositions 13.6 1 and 13.21 we have 



il—ioo 



< Ce- 



9i+ioo 



\X\-'^'il + mL^f\mi.,) ^,MX\ 



LHi) 



< Ce-^'Vlm 



and 



(3.23) 



\s. 



2|l2 



< 



1 

(27r)'^ 


P.V 


^ 1 1 


+ (^ 


!7r)'^ 



-9i+ioo 



— &i— loo 



X-^e^^dX / e^*-«7(a;i,6rfC 



L2 



P.V. 



h+ir 



X-^e^^dX 



i\~ir 



e''<fixuM 



L2 



<Ce-^^*l/lL2, 

by direct computations, noting that the integral in A in the first term is identically zero. 
This completes the proof of the first inequality stated in the proposition. Derivative bounds 
follow similarly. D 
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Remark 3.7. Here, we have used the A^'^ improvement in (j3.5p over (j3.4p together with 
modifications introduced in |KZ| to greatly simphfy the original high-frequency argument 
given in |Z3j for the shock case. 

3.3. Boundary estimates. For the purpose of studying the nonzero boundary perturba- 
tion, we need the following proposition. For h := h{x,t), define 

(3.24) Vh{t):=i\ht\ + \h^\ + \h^^\)it), 
and 

(3.25) Th{t) := I IXY. ^Vk^'^^ + '^^^ ^''' * " *' °' y)^(^' *) ^^^^' 

where G{x, t; y) is the Green function of dt — L. This boundary term will appear when we 
write down the Duhamel formulas for the linearized and nonlinear equations (see p.37p and 
(j4.55p ). Noting that for the outflow case, the fact that G{x,t;0,y) = simplifies Th to 

(3.26) Th{t)= / Gy,{x,t-s;0,y)B^'^hdyds. 

Jo Vm^-i 

Therefore when dealing with the outflow case, instead of putting assumptions on h itself 
as in the inflow case, we make assumptions on B^^h, matching with the hypotheses on 
VF-coordinates. 

Proposition 3.8. Assume that h = h{x,t) satisfies 

\hit)\L,<Eoii+t)-^''+'y\ 

X 

(3.27) \h{t)\Lf<Eo(l + t)-'"^ 

\Vn{t)\^,^^^\,,+-s<E^{l + t)-''l^-\ 

for some positive constant Eq; here \^\ = [{d — l)/2] -|- 2, and e > is arbitrary small for 
d = 2 and e = for d > 3. For the outflow case, we replace these assumptions on h by 
those on B^^h. Then we obtain 

\Th{t)\L2<CEo{l + t)-^''-'^/\ 

(3.28) |r/i(t)|^2,oo < CEoil + t)-('^+l)/^ 



\Thit)\Loo <CEo{l + t) 



'd/2 



and derivative bounds 
(3.29) 



\d^Th{t)\ 2,oo <CEo{l + t) 



-(d+l)/4 



|a|r/i(t)| 2,oo < CEo{i + t)-^'^+^'^/\ 



for all t > 0. 

Proof. We first recall that G{x, t — s; y) is a solution of {ds — Ly)*G* = 0, that is, 

(3.30) _ G. - Y.(GAny, + E G^i. = J2(Gy,B'%- 

3 3 jk 
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Integrating this on R^ x [0, t] against 

(3.31) g{yi,y, s) -.= e~y^h{y, s), 

and integrating by parts twice, we obtain 

^^ = - f I ( E Gy.B''^ + E GA')9y.dyds 
Jo JRi ^y Y ^ 

where, recalhng that 

S{t)f{x)= f Gix,t;y)fiy)dy, 

we get 

- /V E i^y^^'' + E GAAgy^dyds 
Jo Jm^^^^ Y 

= - l"sit-s)[- Y,iB''9.,K +Y.A^g.,)ds 



and 



= - j S{t-s)[gs + Y,^x,9)ds + g{x,t)-S{t)g{x,Q). 
Therefore combining all these estimates yields 
(3.32) Th = g{x,t) -S{t)go- S{t - s){gs - La:g{x,s))ds 



with go{x) := g{x, 0) and L^g = - ^jiA^g)xj + Y.jk{B^''gxk)xj- 

Now we are ready to employ estimates obtained in the previous section on the solution 
operator S{t) = Si{t) + S2{t). Noting that 

\g\a ^ G\h\LV, 
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we estimate 

|r/i|L2 < \9\l^ + \Si{t)go\L2 + \S2{t)go\L2 

+ / \Si{t - s){gs - Lg)\i^2 + \S2{t - s){gs - Lg)\i^2ds 
Jo 

<\h{t)\L2+C{l + t)-^\go\Li+Ce-^'\go\H^ 

X 

+ Al + 1 - s)-^''-^y\\gs\ + \Lg\)Li + e-^(*-'^)(|g,| + {LgD^scis 
Jo 

<\h{t)\L2+C{l + t)-^\ho\LinHi 

JO a: X 

<CEo{l + t)-^ 
and similarly we also obtain 

\Th\^2,^ < \h{t)\L2+C{l + t)-'^\ho\Linm 

X,Xl ^ XX 

(3-33) +C Al + t-s)-('^+i)/^|P.(.)L. +e-^(*-)|P,(.)|^4ds 

JO a: i: 

and 

|r/i|L- < \Ht)\Lf+c{i + t)-l\ho\^,^^i,\+3 

(3-34) +C Al + i-s)-''/'|Ph(5)lLi +e-^(*-^)|P/.(s)l M+3ds 

Jo " ■"« 

Similar bounds hold for derivatives. 

This completes the proof of the proposition. D 

3.4. Duhamel formula. The following integral representation formula expresses the solu- 
tion of the inhomogeneous equation p.lh in terms of the homogeneous solution operator S 
for f,h = 0. 

Lemma 3.9 (Integral formulation). Solutions U of ()3.1I) may be expressed as 

(3.35) Uix, t) =Sit)Uo + [ S{t- s)f{-, s) + TU{0, 5, t) 

Jo 

where U{x,0) = Uo{x), 

(3.36) TU{0,x,t):= f I {y Gy.B^^ + GA^){x,t - s;0,y)U{0,y,s)dyds, 

Jo ilRd-i j 

and G{-,t;y) = S{t)6y{-) is the Green function of dt — L. 
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Proof. Integrating on M^ the linearized equations 

{ds - Ly)U = f 

against G{x, t — s;y) and using the fact that by duality 

ids-LyyG*ix,t-s;y)=0, 
we easily obtain the lemma as in the one-dimensional case (see [YZIINZ] ). recalling that 

S{t)f= [ Gix,t-y)fiy)dy. 



a 

3.5. Proof of linearized stability. 

Proof of Theorem \1.6[ Writing the Duhamel formula for the linearized equations 

(3.37) U{x, t) =S{t)Uo + Th{x, t), 

with Th defined in (|3.25p . where U{x,0) = Uq{x) and U{0,x,t) = h{x,t), and applying 
estimates on low- and high-frequency operators 5i(t) and 52 (t), we obtain 

\U{t)\L2 < \Si{t)Uo\L2 + \S2{t)Uo\L2 + \rh{t)\L2 

(3.38) < C{1 + t)~^\Uo\Li + Ce-'i'lUolH^ + CEo{l + t)-('^-i)/4 
<C{l + ty^{\Uo\Linm+Eo) 

and 

|C/(t)|Loo < \Si{t)Uo\L^ + \S2{t)Uo\L^ + |r/l(t)|ioo 

< C(l + t)-^Uo\Li + C|52(t)C/o|^[(d-i)/2,+2 + CEoil + 1)-''/^ 
(3 39) 

<C{l + t)-2\Uo\L^+Ce-^'\Uo\Hli,-i)m+5+CEo{l + t)-''/^ 

< c(i + ty^ (^\UQ\^i^jj[(d-i)/2]+5 + Eq). 

These prove the bounds as stated in the theorem for p = 2 and p = oo. For 2 < p < oo, we 
use the interpolation inequality between L^ and L°°. D 

4. Nonlinear stability 

4.1. Auxiliary energy estimates. For the analysis of nonlinear stability, we need the 
following energy estimate adapted from |MaZ4l[NZllZ4] . Define the nonlinear perturbation 
variables U = (n, v) by 

(4.1) U{x,t):=U{x,t)-U{xi). 

Proposition 4.1. Under the hypotheses of Theorem \1.'T\ let Uq G H^ and U = {u,v)'^ be a 
solution of (jl.2p and ()4.ip . Suppose that, for < t <T, the Wx'°^ norm of the solution U 
remains bounded by a sufficiently small constant C > 0. Then 

(4.2) \U{t)\jjs < Ce-''\Uo\h + cJ\-'^'-^\\U{t)\1, + \B,ir)\')dT 
for all < t < T, where the boundary term Bh is defined as in Theorem \1.7\ 
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Proof. Observe that a straightforward calculation shows that |C/|h'- ~ |VF|_h''-, 

(4.3) W = W -W ■.= W{U)-W{U), 

for < r < s, provided |C/|i4/2,oo remains bounded, hence it is sufficient to prove a corre- 
sponding bound in the special variable W. We first carry out a complete proof in the more 
straightforward case with conditions (Al)-(A3) replaced by the following global versions, 
indicating afterward by a few brief remarks the changes needed to carry out the proof in 
the general case. 

(Al') A^W),A°,Al^ are symmetric, i° > 6*0 > 0, 

(A2') Same as (A2), 

(A3') I^ = (|/,) , &'' = B'^^ = (^^ ~^,y EOe.^'>^leP,andG^0. 
Substituting (j4.3p into (II. 4p . we obtain the quasilinear perturbation equation 

(4.4) A°Wt + J^ Am^. = Y,{B'''W^,)x, + M^W^, + ^(M^T^, J,, 

J jk j 

where A^ := A^{W + W) is symmetric positive definite, A^ := A^(W + W) are symmetric, 
Ml = A\W + W)- A^{W) = ( I dA\W + eW)d9)w, 

As shown in |MaZ4j . we have bounds 

(4.5) \A^\ < C, |A0| < C\Wt\ < C{\W^\ + \wii\) < CC, 

2 

(4.6) \d,A'\ + \dlA'\ < C{J2\d'.W\ + \W,,\) < CiC+\W,,\). 

k=l 

We have the same bounds for A^ , B^ , and also due to the form of Mi , M2 , 

(4.7) \Mi\,\M2\<C{C + \W^A)\W\. 

Note that thanks to Lemma 11.31 we have the bound on the profile: \Wxj^\ < Ce~^'^^', as 
xi — > +00. 

The following results assert that hyperbolic effects can compensate for degenerate vis- 
cosity B, as revealed by the existence of a compensating matrix K. 

Lemma 4.2 ( [KShj ). Assuming (Al'), condition (A2') is equivalent to the following : 

(Kl) There exist smooth skew- symmetric "compensating matrices" K{^), homogeneous 
degree one in ^, such that 

(4.8) ^e^Y.^j^kB^'^ - KiO{A')-'Y,^kA'')iW+) > ^sICP > 

j,k k 

for all C £ R'^ \ {0} . 
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Define a by the ODE 

(4.9) a^, = -sign{A\^)c,\W^,\a, a(0) = 1 
where c* > is a large constant to be chosen later. Note that we have 

(4.10) {a^Ja)A\^ < -c,0i|Ty^J =: -uj{xi) 
and 

(4.11) \a^Ja\ < c^lVF^J = e^^uj{xi). 

In what follows, we shall use (•, •) as the a- weighted L^ inner product defined as 

{f,g) = (a/,5)L2(Kd) 
and 

4=0 \a\=i 

as the norm in weighted H^ space. Note that for any symmetric operator S, 
{SU,,f) = -^{S,J,f), i/1 
{Sf.,J) = -Ins., + {a.Ja)S)f,f) - ^{Sf,f)„ 

where (•, •)g denotes the integration on Mg := {xi = 0} x M . Also we define 

\\f\\o,s = ||/|lH=(Mg) = E E {d^f^d^iy^'- 

i=0 \a\=i 

Note that in what follows, we shall pay attention to keeping track of c*. For constants 
independent of c*, we simply write them as C. Also, for simplicity, the sum symbol will 
sometimes be dropped where it is no confusion. We write ||/a;|| = ^- \\fxj\\ and ||i9^/|| = 

Yl\a\=k W^xfW- 

4.1.1. Zeroth order "Friedrichs-type" estimate. First, by integration by parts and estimates 
K5\\ . (foil , and then KW\\ . we obtain for j ^ 1, 

-{Amx^,W) = ^{AiW,W)<C{iC + \Wx,\)w',w')+C\\w''\\l 



and for j = 1, 

^{{Al^ + {axJa)A')W,W) + ^ 

1 

2 
1 



{A'Wx,,W) = -{{Al^ + (axJa)A')W,W) + -{A'W,W), 

< ^{iaxJa)AW,w') + C((C + \Wx,\)\W\ + u;ixi)\w''\,\W\) + J^ 



< -l{u;{x)w',w') + cue + \Wx,\)w',w') + C{c.)\\w''\\l + Jl 
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where J^ denotes the boundary term ^{A^W,W)q. The term {\Wxi\w^ ,w^) may be easily 
absorbed into the first term of the right-hand side, since for c^, sufficiently large, 

(4.12) {\W,,,\w',w') < {c,eir\Lo{xi)w',w') < -^{io{xi)w',w'). 

Also, integration by parts yields 

J/||2 , r^L ,/_ \„..n „„II\ lulk„,TI „„II\ 






where we used the fact that B^^Wx-W = IP^wl} -w^^ , noting that B has block-diagonal form 
with the first block identical to zero. Similarly, recalling that M2 = B^^iW + W) — B^'^iW), 
we have 

< C{\WxA\W\,\w'J\) + C{lo{xi)\W\,w'') - {mlWx„w'')^ 

< eikf llo + C(^e{u;ix,)w',w') + C ic,)\\w' ' \\l) - {mlW^^w''), 

for any small ^, e. Note that C is independent of c^,. Therefore, for ^ = 9/2 and c* sufficiently 

large, combining all above estimates, we obtain 

(4.13) 

— {A'^W, W) = {A^Wt, W) + ^{A^tW, W) 

= {-A^x^ + {B^'^WxJx, + MiWx, + {MiWx,)x,,W) + i(^?VF, W) 
< -h^{u:{xr)w' ,w') + e\\wi%] + CC\W\\l + C{c,)\W'\\l + 1° 
where the boundary term 

(4.14) It := 1{A'W,W), - {b''wil,w'')^ - {mlWx„w'% 

which, in the outflow case (thanks to the negative definiteness of ^j;;^), is estimated as 

(4.15) 4° < -||k^||g,o + CiWWlo + \\wi'\\o,oW\\o,o), 
and similarly in the inflow case, estimated as 

(4.16) 4° < Ci\\W\\lo + \\w'J \\o,oW\\o,o). 
Here we recall that || . ||o,s := || • H/fsmd-j. 
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4.1.2. First order "Friedrichs-type" estimate. Similarly as above, we need the following key 
estimate, computing by the use of integration by parts, (14.12p . and c=k being sufficiently 
large, 

j 

j 

We deal with the boundary term later. Now let us compute 

(4.18) ^±{A^w,^,W,^) = {W,^,iA'WtU) - {W,^,AlWt) + ^{A^,W,^,W,^). 

We control each term in turn. By ()4.5p and ()4.6p . we first have 

and by multiplying {A^)~^ into (14. 4p . 

\{W,^,AlWt)\<C{{C + \W,AW.li\WA + \wii\ + \W\)) 

<^\\wii\\l + C{{C + \W,,\)wi,wi) + C{{C + \W,,\)w\w')+C\\w''\\l 

where the term {\Wxi\wl,,wl.) may be treated in the same way as was {\Wxi\w^ ,w^) in 
()4.12p . Using (j4.4p . we write the first term in the right-hand side of (j4.18p as 

{Wx,,{A^Wt)x^ ={Wx^, [-A^x^ + (B^'^Wx,).^ + MiWx, + (M^'VF^J^J^J 
= - {W,,,Am,^,^) + {W,^,-AiW.^ + (M.W,,),^) 

- {{a^Ja)W^^, [{B'^W,,).^ + (MJW^.J.J) 



< 



1 

4L 

C 



{uj{xi)wi,wi) + 9\\wii\\l 



+n 



where /^ denotes the boundary terms 



(4.19) 



/I : = ^{W^^,A'W^^), - {W,^, [(B^^T^.,)., + (M^W^.U), 



1 



^{W,^,A'W,^), - «, [(6^^<)., + (m^iy.J^J)^, 

and we have used (A3) for each fixed i and S,j = {Wxi)x- to get 

(4.20) 5^(H^x..,,i?^''H^x,.,) >0Y,\\W,,.,\\l 

jk j 
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and estimates ()4.17p . (|4.12p for w^ ,w^^, and Young's inequality to obtain: 

{W^, -AiW, + (MiVF,.J,) < C((C + \w^AWc.\, \w^\ + \w\). 

-ewwiiwl + c{\wii\ + ujixi)\w'Jiic + l^xiDkf I) 

C{\wii\ + u;{x,)\w'J\,{C + \W^A)i\WA + \W\)). 
Putting these estimates together into (|4.18p . we have obtained 
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^±{A'w,,w,) + \e\\wii\\l + \Hx,)wi,wi) 



(4.21) 



< c 



I „..I\ 



,//||2 



C\\w'\\i + {\W,,\w',w') + C{c,)\\w''\\i 



+Il 



Let us now treat the boundary term. First observe that using the parabohc equations, 
noting that A^ is the diagonal-block form, we can estimate 

(4.22) iV'wi%^iO,x,t)<c{\wi'\ + \W,^\ + \W\yO,x,t) 

and thus for i ^ 1 

< 



(4.23) 



\wi^x\i\W\ + \wii\ 



< c 



P^P + |y;f |2 + |u;|||2 



and for i = 1, using b = 6 , ()4.22p . and recalling here that we always use the sum 
convention, 



E(^''0-i = liib^'^u + (b^'O., + bllwii - ^< 



(4.24) 



Therefore 






<C{\wi^\ + \W\ + \W^^\ + \w:, 



{wil,[{b''wilU + {mlW^,U])o 



< e 



\wi\^ + c 



.JI\2 



Wt 



+ \W\' + \wi'\' + \wii 



^xx\ 



For the first term in /^, we consider each inflow/outflow case separately. For the outflow 
case, since ^4}^ < —9i < 0, we get 

A^W..W.<-^-^\wif + C\wi'\'. 
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Therefore 

(A oc:\ rl ^ ^1 /" I ^|2 I /" fljxrl2 I I //|2 , i II\2 , i II\2 

Meanwhile, for the inflow case, since A}y^ > ^i > 0, we have 

In this case, the invertibility of A}y^ allows us to use the hyperbolic equation to derive 

\wi^\<Ci\w\\^\wi'\^\wi\). 
Therefore we get 

(4.26) ll<\ (|T^|2 + |iy,|2 + |^||2 + |^//|2 + |^//|2\ 

Now apply the standard Sobolev inequality 

(4-27) k(0)P < C||w||j^2(K)(||u;3;||i2(R) + ||'u;||l2(r)) 

to control the term |tL'^{(0)p in l\ in both cases. We get 

(4.28) / \wi{^<^\\wii\g^C\\wi'\g. 
Using this with e' = 61/8, (liHI) . and ([425]), the estimate (1421]) reads 

(4.29) «* 

where the (new) boundary term /^ satisfies 

for the outflow case, and 

(4.31) ll < [ (\W\^ + \Wt\^ + \W^\^ + \wg\^ 

for the inflow case. 



4.1.3. Higher order "Friedrichs-type" estimate. For any fixed multi-index a = (a^^ , • ' ' > ckx^)) 
ai = 0, 1, \a\ = k = 2, ...,s, by computing ■^{A'^d^Wjd^W) and following the same spirit 
as the above subsection, we easily obtain 

(4.32) fc-i 

i=l 
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where 

aa . _ D«i . . . fjOd f)a+l — \^ ffOi . . . flf^df) f)i — Sr^ fi0l . . . f)f^d 

^x ■ — ^xi ^Xd> ^x ■— Z^^xi ^xa'-'xj: ^x — 2.^ '~'x\ ^x^ 

J \l3\=i 

and the boundary term /" satisfies 

(4.33) /r<-y/ \9>'\' + C ( J] |9>^^|^ + ^ |5>^P + ^ |5|..^^|^ 
for the outflow case, and 

. k [{k+l)/2] k 

(4.34) I^< l[^\dy\'+ Y. \diw''\' + Y\dlW\ 

"^*0 1=0 i=l i=0 

for the inflow case. 

Now for a with ai = 2, ..., s we observe that the estimate (|4.32p still holds. Indeed, using 
integration by parts and computing j^{A^d^W,d^W) as above leaves the boundary terms 
as 

(4.35) I^ := ^{d^W,A'd^W), - {d^w^^ ,d^[ib^'^wii) + (mlW,,)])^. 
Then we can use the parabolic equations to solve 

Using this we can reduce the order of derivative with respect to xi in d" to one, with the 
same spirit as ()4.23p and ()4.24p . Finally we use the Sobolev embedding similar to (j4.28p to 
obtain the estimate for the normal derivative dxi , and get the estimate for I^ as claimed in 
([03]) and gJH). 

We recall next the following Kawashima-type estimate, presented in |Z3] . to bound the 
term Htf^H^, appearing on the left hand side of (j4.32p . 

4.1.4. "Kawashima-type" estimate. Let K(£) be the skew-svmmetrv in (j4.8p . Using Plancherel's 
identity and the equations (j4.4p . we compute 



(4.36) 



where 



(4.37) 



■■{iK{om'w,iiCYWt) 

+ {iK{0{iCYW,{iCrH), 



H:=Y, («)-'4 - (A'^r'AAw, 



j 



+ iA''r\Y{B^'W,Jx,+M,W,,+Y,iMiWx,] 

jk j 
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By using the fact that |(^^)"M^ - (A'^y^A^ = 0{C + \W^i\), we can easily obtain 

r+l 
k=0 

Meanwhile, applying (j4.8p into the first term of the last line in ()4.36p . we get 

myw, - K{0{Al)-'^C,Aim''W) 

j 

>omr'w\\i-cmr'w''\\i 

Putting these estimates together into ()4.36p . we have obtained the high order "Kawashima- 
type" estimate: 

^^{Kid,)diw, diw) < - e\\di+W\\l + c\\w''\\^.^2 

(4.38) r+l 



4.1.5. Final estimates. We are ready to conclude our result. First combining the estimate 
(f429]) with (f4T3]) . we easily obtain 

^^^(^{A'^W^W,) + M{A'W,W) 



< - [-^WwiiWl + -{u;{xi)wi,wi)j 

+ C[c\\w'\\l + {\W,,\w',w') + C{c,)\\w''\\l) + 4^ 



- ^(^{u;{xi)w', w') + 9\\wi'\\l) + CMCWw'Wl + MC{c.)\\w''\\l + Mlj^ 

By choosing M sufficiently large such that M6 » C(c*), and noting that c*0i|Wa;J < 
tL'(xi), we get 

i^ ((AOt^,, W^) + M{A^W, W) 

(4.39) < - (^e\\w''\\l + {u;ixi)w',w') + {u;ixi)wi,wi)) 

+ C[C\W\\1 + Cic,)\W'\\l) + n + Mil 

We shall treat the boundary terms later. Now we use the estimate (j4.38p (for r = 0) to 
absorb the term HSxif^^Ho into the left hand side. Indeed, fixing c* large as above, adding 
(j4.39p with (j4.38p times e, and choosing e, C sufficiently small such that eC(c*) <C 0, e <C 1 
and C <^ 602, we obtain 
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^^ ((^°H^,., W^) + M(^°M^, W) + e{KW,, W) 



+ c[c\W\\i + Cic.)\\w^^\\l)-^-fU\\l 

+ Ce(^\\w''\\l + CWw'Wl + {u;{xi)w',w') + {u;{xi)wi,wi)) + ll + MI^, 



2 

where h ■= ll +MI^. 

In view of boundary terms I^ and /^, we treat the term /^ in each inflow/outflow case 
separately. Recah the inequahty (j4.28p . ||w^.{||o,o ^ C\\w^^\\2- Thus, using this, for the 
inflow case we have 

(4.40) /o < C{\\W\\l, + ||^,^^||o,olk''||o,o) < C{\\W\\l, + \W^\\l, + 6||u;^^i) 
and for the outflow case, 

(4.41) ^ 

S — ir\\w ||o,o + '-"UI^ llo,o + ll^x llo,o + ^ll^ \\2)- 

Therefore these together with ()4.30p and (j4.3ip . using the good estimate of ||tu^^||o, yield 

(4.42) h<-%l iW? + \wi\') + C [ (\w^^\' + l^f |2 + \wi^\' + \wg^' 

for the outflow case, and 

(4.43) ll<l (|W^P + |Wt|2 + |W^s|2 + |^i;f |2' 



for the inflow case. 

Now by Cauchy-Schwarz's inequality, |i^(^)| < C|^|, and positive definiteness of A^, it is 
easy to see that 

(4.44) £ : = (A'^W^^W^) + M{A°W,W) + e{K{d^^)W,W) ~ ||VF||^i ~ llW^ll^i- 

The last equivalence is due to the fact that a is bounded above and below away from zero. 
Thus the above yields 

piwm < -03£{w){t) + c{c.)[\\wmh + \si{t)f), 

for some positive constant ^3, which by the Gronwall inequality implies 

(4.45) ||W^(t)|||,. <Ce-^*||T^o||?,i+C(c.)^*e-^(*-^)(ir(r)||i. + |^i(r)|2)dr, 
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where W{x,0) = Wo{x) and 

(4.46) \Bi{T)\^:=f (\W\^ + \Wt\^ + \Wi\^ + \wgf 

for the inflow case, and 

(4.47) \B,{t)\':=[ f|«;^^|2 + |^|^|2 + |^//|2 + |^//,2 



-^xx\ 



for the outflow case. 

Similarly, by induction, we can derive the same estimates for W in H^ . To do that, let 
us define 

£i{W) := {A°W^, W^) + M{A°W, W) + e{KW^, W) 

£k{W) := {A^d^W,d^W) + M£k-i{W) + e{Kd^W,d^-^W), k < s. 

Then similarly by the Cauchy-Schwarz inequality, £s{W) ~ ||W^|||/s) and by induction, 
we obtain 

ps{Wm < -e,£,{wm + C{cMW{t)\\l, + \Bf,{t)\% 

for some positive constant ^3, which by the Gronwall inequality yields 



|2 / ry„-et 



\Wo\\h+Cic.) f e-'^'-^H\\W{T)\\l, + \B,{r)\')dT, 
Jo 



(4.48) \\W{t)rH.<Ce^ 

'0 

where W{x,0) = Wo{x), and Bh are defined as in ()1.14p and (J1.15p . 

4.1.6. The general case. Following |MaZ4t [Z3] . the general case that hypotheses (Al)-(A3) 
hold can easily be covered via following simple observations. First, we may express matrix 
A in (|ill|) as 

(4.49) A^(^W + W)=A^ + iC + \W,,\)(^JI^^^ ^j}i 

where A^ is a symmetric matrix obeying the same derivative bounds as described for A^ , 
A^ identical to A^ in the 11 block and obtained in other blocks kl by 

AliiW + W) = AliW) + AliW + W)- AliW) 

(4.50) =Ali{W+) + 0{\W,\ + \W,,\) 

and meanwhile, A^ , j ^ I, obtained by A^ = A^{W+) + 0{C + \Wxi\), similarly as in (|4.50p . 
Replacing A^ by A^ in the k^^ order Friedrichs-type bounds above, we find that the 
resulting error terms may be expressed as 

{d',o{c + \w,,\Wl\d',+W\), 

plus lower order terms, easily absorbed using Young's inequality, and boundary terms 

k 

o(j;ia>^^(o)pV(o)|) 

i=0 
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resulting from the use of integration by parts as we deal with the 12-block. However these 
boundary terms were already treated somewhere as before. Hence we can recover the same 
Friedrichs-type estimates obtained above. Thus we may relax (^1') to (Al). 

Next, to relax (^3') to (^3), first we show that the symmetry condition B^'^ = B^^ is 
not necessary. Indeed, by writing 

jk jk jk 

we can just replace B^^ by B^^ := ^{B^^ + B^^), satisfying the same {A3'), and thus still 
obtain the energy estimates as before, with a harmless error term (last term in the above 
identity). Next notice that the term g(Wx) — giWx^) in the perturbation equation may be 
Taylor expanded as 

MWx,Wx,)+gi{Wx„Wx)J ^ \0{\Wx\- 



dlw,dl , ,..,.., ,2- 



The first term, since it vanishes in the first component and since \Wx\ decays at plus spatial 
infinity, yields by Young's inequality the estimate 

which can be treated in the Friedrichs-type estimates. The (0,O(|W2:p) nonlinear term 
may be treated as other source terms in the energy estimates. Specifically, the worst-case 
term 



P{Wx\ 

= -{dl-^W ,dl-^0{\Wx\^)) - dW\^)dl-'0{\Wx\^m 
may be bounded by 

\\dx'^'^w^'\\L2\\W\\w2,^\\W\\Hk - d';w^\'d)d^x~^0{\Wx\'^){^). 

The boundary term will contribute to energy estimates in the form (j4.35p of /", and thus 
we may use the parabolic equations to get rid of this term as we did in (j4.23p . (j4.24p . Thus, 
we may relax {A2>') to (A3), completing the proof of the general case (^41) — (A3) and the 
proposition. D 

4.2. Proof of nonlinear stability. Defining the perturbation variable U := U — U , we 
obtain the nonlinear perturbation equations 

(4.51) Ut-LU = Y,QHU,Ux)x,, 

j 
where 

Q^{U,Ux) = 0{\U\\Ux\ + \U\'') 

(4.52) Q'{U,Ux)x,=0{\U\\Ux\ + \U\\Uxx\ + \Ux?) 

Q'iU,Ux)x,x, = 0{\U\\Uxx\ + \Ux\\Uxx\ + \Ux\^ + \U\\Uxx.\) 
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SO long as |[/| remains bounded. 

For boundary conditions written in [/-coordinates, (B) gives 



h = h-h = {W{U + U)- T^(C/))(0, X, t) 

= {dw/diJ){Uo)u{<d,x,t) + o{\u{<d,x,t)\^). 

in inflow case, where {dW / dU){UQ) is constant and invertible, and 

h = h-h = {w^\U + U)- w^^{U)){0, X, t) 

= {dw''/dtj){U^)U{^,i,t) + 0{\U{^,±,t)\'') 
(4.54) 

= m{hi 62)(C/o)^(0,x,t) + 0(|t/(0,x,t)|2) 

= m5(i7o)[/(0, X, t) + 0(1^/(0, X, t)p) 



for some invertible constant matrix m. 
Applying Lemma 13.91 to (j4.5ip , we obtain 

(4.55) Uix,t)=S{t)Uo+ f S{t-s)Y,d.,Q'{U,U,)ds + TU{Q,x,t) 

where ^(x,0) = Uq{x), 

(4.56) TU{0,S:,t):= f [ (TGyW^ + GA^){x,t-s;Q,y)U{Q,y,s)dyds, 

Jo JKd-i ^ 

and G is the Green function oi dt — L. 
Proof of Theorem |j.7| Define 

at) := supf |[/(s)|i2 (1 + s)"-^ + \U{s)\l^{1 + s)i 



(4.57) 

+ {\U{s)\ + \U,{s)\ + |5|C/(s)|)^2,^ {l + s)^ 

We shall prove here that for all t > for which a solution exists with C(t) uniformly 
bounded by some fixed, sufficiently small constant, there holds 

(4.58) m < C{\Uo\L^nH^ + Eo + C{t)^). 
This bound together with continuity of C{t) implies that 

(4.59) C{t)<2G{\Uo\L^nH^+Eo) 

for t > 0, provided that |f/o|LinH= + Eq < 1/4C^. This would complete the proof of the 
bounds as claimed in the theorem, and thus give the main theorem. 

By standard short-time theory/local well-posedness in H^, and the standard principle 
of continuation, there exists a solution U E EP^ on the open time-interval for which |C/|h'' 
remains bounded, and on this interval ({t) is well-defined and continuous. Now, let [0,T) 
be the maximal interval on which |C/|hs remains strictly bounded by some fixed, sufficiently 
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small constant 6 > 0. By Proposition l4.lt and the Sobolev embeding inequality |C/|i4/2,oo < 
CIUIh", we have 



(4.60) 



\Uit)\h < Ce-''\Uo\ls+C j\-'^'-^^(\U{T)\l, + \Bn{Tt)dT 



<C{\U^\ls+El + atf){l + t)- 



-{d-l)/2 



and so the solution continues so long as C, remains small, with bound ()4.59p . yielding exis- 
tence and the claimed bounds. 

Thus, it remains to prove the claim (j4.58p . First by ()4.55p . we obtain 



(4.61) 



\U{t)\L2 <|5(t)f/o|i2 + / |5i(t - s)d^^Q^{s)\L2ds 

Jo 

+ j \S2{t - s)d^^Q^ {s)\L2ds +\TU{'d,X,t)\L2 

</i + /2 + i3 + |rc/(o,x,t)|i2 



where 



h : = |5(t)C/o|z.2 < C(l + t)-^|C/o|Lin//3, 
h-= I \Si{t - s)d^,^Q\s)\L2ds 

<C (l + t-s)'V-5|Qi(s)|ii + (l + s)-^|QJ'(s)Li,oods 

JO "'•"'1 

<C {l + t-sy'^-'^\U\]j^ + {l+t-s)-^(\U\]2,o. +|C/x|%,oo ]ds 



'Li'' 



X ,X^ 



<Ci\Uo\h+C{t)' 



[l+t- S) 4 2 



(1 + s)- — 



+ (l + t-s)"^(l + s)"^ 



ds 



<cii+t)-^{\uo\is + atf) 
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and 

^3:= / \S2{t-s)d^^Q^{s)\L2ds 
t 



Jo 
Jo 

Jo 

< C{\Uo\ls + C{t?) f e-^(*-^)(l + sy'^ds 
Jo 



<c{i+t)-—{\uo\is+atf). 

Meanwhile, for the boundary term \TU{d,x,t)\i2, we treat two cases separately. First for 
the inflow case, then by (j4.53p we have 

|C/(0,x,t)|<C|Mx,t)|+O(|C/(0,x,t)|2), 

and thus |[/(0,x,t)| < C\h{x,t)\, provided that \h\ is sufficiently small. Therefore under the 
hypotheses on h in Theorem II. 7t Proposition 13.81 vields 

\TU{Q,;-)\Ll<CE^{l+t)-'-^. 

Now for the outflow case, recall that in this case G{x, t; 0, y) = 0. Thus ()4.56p simplifies 
to 

(4.62) TU{0,x,t)= / Gy^{x,t-s;0,y)B^^U{0,y,s)dyds. 

Jo jRd-l 

To deal with this term, we shall use Proposition 13.81 as in the inflow case. In view of (|4.54p . 

\B^'U{0,y,s)\<C\h{y,t)\+O{\U{0,y,s)\^), 

and assumptions on h are imposed as in Theorem 11.61 so that (j3.27p is satisfied. To check 
the last term O{\U{0)\'^), using the definition (j4.57p of C{t), we have 

\O{\U{0,y,s)\^)\L2 <C\U\l^\U\^2,.. <CC'(t)(l + s)-i-^ 

\O{\U{0,y,s)\^)\L^ < C\U\l^ < Cemi + s)-'' 

and for the term Vh with h replaced by O{\U{0, y, s)p), using the standard Holder inequality 
to get 

|2?,.|^[,.-l)/2, + 5 < C|[/U.o|[7|^. < CC'(t)(l + 5)-'^/2-('^-l)/4. 

X 

We remark here that Sobolev bounds (|4.60p are not good enough for estimates of Vh in L^, 
requiring a decay at rate (1 + t)~'^''^~^ for the two-dimensional case (see Proposition 13. 8p . 
This is exactly why we have to keep track of Uxx in i>^'°° norm in (^{t) as well, to gain a 
bound as above for Vh. 
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Therefore applying Proposition [3]8l we also obtain (j4.62p for the outflow case. Combining 
these above estimates yields 

(4.63) \u{t)\L2il+t)'^ <C{\Uo\L^nHs+Eo + Citf)- 

Next, we estimate 

|C/(i)|^2,oo <|5(t)C/o 1^2.00 + / |5i(t-s)a^.Q^'(s)|^2,oods 

x,x^ x.x^^ I r\ x,x-^ 

^^■^^) + / |cS2(t-s)5^^Q^(s)|^2,oods + |rC/(0,x,i)|^2.oo 

< Ji + J2 + J3 + |rc/(o, X, t)L2,oc 

where 

Ji : = \S{t)Uo\j2,^ < C(l + t)-'^|f/o|Lini?4 

J2 ■ = \Si{t - s)dxjQ^{s)\j^2,oa ds 

Jo ^'^1 

Jo ^'^1 

^0 X,X^ XjX-^^ / 

JO 

d+1 ^ ^ d+1 

+ (l + t-s) ^(1 + s) 2- (is 

<C(l + t)~^(|t/o|?,.+C(t)') 
and (by Moser's inequality) 

J3 ■= \<S2{t - s)5a;.Q-'(s) I ^2,00 ds 

Jo ^'^1 

<C [ e-'^'~'^\d,^Q^is)\H,ds 
Jo 

<C [\-'^'~'^\U\L^\U\H^ds 

Jo 

<C{\Uo\h+Citf) / e"'^*"^)(l + s)-f(l + s)-Vds 
Jo 

KCil + ty'i^QUolh+atf)- 

These estimates together with similar treatment for the boundary term yield 



d+l 



(4.65) |C/(t)l r?,oo (l + t)— < Ci\Uo\LinH^ + Eo + ((t) 



2n 



X ,X^ 
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Similarly, we have the same estimate for \Ux{t)\j 2,00 . Indeed, we have 



(4.66) 



\Ux{t)\r'2,oo <|9x5(t)C/o|r2,oo + / \dxSi{t- s)dx.Q\s)\r'2,<^ds 

x,x-^ x,x-^ I r\ x^x-^ 

If) X,Xl^ X,3. 

<Ki+K2 + K3 + \d,TUiO,x,t)\,2, 



x,x-i 



where K2 and K^ are treated exactly in the same way as the treatment of J2 , <^3 ; yet in 
the first term of K2 it is a bit better by a factor t~^'^. Similar bounds hold for \Uxx\ in 
L^'°°, noting that there are no higher derivatives in xi involved and thus similar to those 
in (fi^Hi]) . 

Finally, we estimate the L°° norm of U . By Duhamel's formula (j4.55p . we obtain 



\U{t)\Loo <\S{t)Uo\Lo^ + j \Si{t-s)dx^Q^{s)\Lo^ds 



(4.67) 



+ / \S2{t-s)dx^Q^{s)\L^ds + \TU{{),x,t)\L^ 
Jo 

<Li + L2 + L3 + \TU{0,x,t)\L^ 



where the boundary term is treated in the same way as above, and for I7I = [{d— l)/2] + 2, 



Li : = \S{t)Uo\L^ < C(l + t) 2 |t/o|iin^bi+3, 

L2- = \Si{t - s)dx^Q\s)\Loods 

ft 



<C I {1 + t - s)-i--2\Q^ {s)\li + {1 + s)~i\Q^ {s)\ri,^ ds 

Jo ^.^1 

<C I {1 + t - syi-^lUljj^ + {1 + t - s)-U\U\'^.2.o. +\Ux\'i 2,0.) ds 



<Ci\Uo\h+Citf) 



t r- 



[l+t-s)-^-^(l + s) 



+ (l + t-s)~2(i + s)- 
<Cil + tr-2{\Uo\h+Citf) 



d+l 



ds 
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and (again by Moser's inequality), 

L3-= 1^2 (t - s)dx^Q^{s)\Lo^ds 

< \S2{t- s)dxjQ\s)\jj\^\ds 

Jo 

<C [ e-^(*-^)|C/|Loo|C/|^|,|+5ds 
Jo 

< Ci\Uo\h + C{tf) f e-^(*-)(l + .)-i(l + sr'-^ds 

Jo 

<C{l + t)-ii\Uo\ls+at?). 
Therefore we have obtained 

(4.68) \U{t)\L^{l + t)i<C{\U^\L^^H^+E^ + at?) 

and thus completed the proof of claim (j4.58p . and the main theorem. D 



Appendix A. Physical discussion in the isentropic case 

In this appendix, we revisit in slightly more detail the drag-reduction problem sketched in 
Examples 1 1 . lf[L2l in the simplified context of the two-dimensional isentropic case. Following 
the notation of [GMWZ5J . consider the two-dimensional isentropic compressible Navier- 
Stokes equations 

(A.l) pt + (pu)^ + {pv)y = 0, 

(A.2) {pu)t + {pu^)x + {puv)y +px = {2p + r])uxx + puyy + {p + r])vxy, 

(A.3) {pv)t + {puv)x + {pv'^)y +Py = fJ-Vxx + (2/J. + 1])Vyy + (p + r])Uyx 

on the half-space y > 0, where p is density, u and v are velocities in x and y directions, 
and p = p{p) is pressure, and p > \r]\ > are coefficients of first ("dynamic") and second 
viscosity, making the standard monotone pressure assumption p'{p) > 0. 

We imagine a porous airfoil lying along the x-axis, with constant imposed normal velocity 
'u(O) = V and zero transverse relative velocity u{0) = imposed at the airfoil surface, and 
seek a laminar boundary-layer flow {p,u,v){y) with transverse relative velocity Uoo a short 
distance away the airfoil, with \V\ much less than the sound speed Cqo and \uoo\ of an order 
roughly comparable to Cqo- 

A.l. Existence. The possible boundary-layer solutions have been completely categorized 
in this case in Section 5.1 of [GMWZS] . We here cite the relevant conclusions, referring to 
[GMWZS] for the (straightforward) justifying computations. 
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A. 1.1. Outflow case (V < 0). In the outflow case, the scenario described above corresponds 
to case (5.15) of jGMWZS] . in which it is found that the only solutions are purely transverse 
flows 

(A.4) {p,v)^{po,V), n(y)=t.oo(l-e^«^^/^), 

varying only in the tranverse velocity u. The drag force per unit length at the airfoil, by 
Newton's law of viscosity, is 

(A.5) tJ'Uy\y=0 = UooPoo\V\, 

since momentum m := p^V = PcxiV is constant throughout the layer, so that {poo, Uoo being 
imposed by ambient conditions away from the wing) drag is proportional to the speed \V\ 
of the imposed normal velocity. 

A. 1.2. Inflow case (V > 0). Consulting again |GMWZ5j (p. 61), we find ioi V > with 
specified {p,u,v){0) of the orders described above, the only solutions are purely normal 
flows, 

(A.6) n = n(0), {p,v) = ip,v)iy), 

varying only in the normal velocity v. Thus, it is not possible to reconcile the velocity u{0) 
at the airfoil with the velocity Uoo >> c some distance away. 

As discussed in |MN| . the expected behavior in such a case consists rather of a combi- 
nation of a boundary-layer at y = and one or more elementary planar shock, rarefaction, 
or contact waves moving away from y = 0: in this case a shear wave moving with normal 
fluid velocity V into the half-space, across which the transverse velocity changes from zero 
to Uoo- That is, a characteristic layer analogous to the solid-boundary case detaches from 
the airfoil and travels outward into the flow field. In this case, one would not expect drag 
reduction compared to the solid-boundary case, but rather some increase. 

A. 2. Stability. If we consider one-dimensional stability, or stability with respect to per- 
turbations depending only on y, we find that the linearized eigenvalue equations decouple 
into the constant-coefficient linearized eigenvalue equations for (p, v) about a constant layer 
(p, v) = {po, V), and the scalar linearized eigenvalue equation 

(A. 7) Xpu + muy = puyy 

associated with the constant-coefficient convection-diffusion equation put + mu = puyy, 
m := pv = poV, p = Po- As the constant layer {pQ,V) is stable by Corollary 11.51 or direct 
calculation (Fourier transform), and ()A.7p is stable by direct calculation, we may thus 
conclude that purely transverse layers are one-dimensionally stable. 

Considered with respect to general perturbations, the equations do not decouple, nor do 
they reduce to constant-coefficient form, but to a second order system whose coefficients are 
quadratic polynomials in e^''^. It would be very interesting to try to resolve the question 
of spectral stability by direct solution using this special form, or, alternatively, to perform 
a numerical study as done in |HLyZ2| for the multi-dimensional shock wave case. 

Remark A.l. For general laminar boundary layers {p,u,v){y), the one-dimensional sta- 
bility problem, now variable-coefficient, does not completely decouple, but has triangular 
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form, breaking into a system in {p,v) alone and an equation in u forced by {p,v). Stabil- 
ity with respect to general perturbations, therefore, is equivalent to stability with respect 
to perturbations of form {p,0,v) or (0,m, 0). For perturbations {p,u,v) = (0,^, 0), the u 
equation again becomes ()A.7p . with p, m still constant, but p varying in y. Taking the real 
part of the complex L^ inner product of u against (|A.7p gives 

3f?A||n|||2 + ||%|li2 = 0, 

hence for ?R.X > 0, n = constant = 0. Thus, the layer is one-dimensionally stable if and only 
if the normal part {p,v) is stable with respect to perturbations {p,v). Stability of normal 
layers was studied in [CHNZj for a 7-law gas p{p) = ap'^ , 1 < 7 < 3, with the conclusion 
that all layers are one-dimensionally stable, independent of amplitude, in the general inflow 
and compressive outflow cases. Hence, we can make the same conclusion for full layers 
{p,u,v). In the present context, this includes all cases except for suction with supersonic 
velocity \V\ > Cqo, which in the notation of |CHNZ| is of expansive outflow type (expected 
also to be stable, but not considered in [CH NZ]). since \v\ is decreasing with y, so that 
density p (since m = pv = constant) is increasing. 

A. 3. Discussion. Note that we do not achieve by subsonic boundary suction an exact 
laminar flow connecting the values {u,v) = {0,V) at the wing to the values (woojO) of the 
ambient flow at infinity, but rather to an intermediate value {uoo,V). That is, we trade 
a large variation u^o in shear for a possibly small variation V in normal velocity, which 
appears now as a boundary condition for the outer, approximately Euler flow away from 
the boundary layer. Whether the full solution is stable appears to be a question concerning 
also nonstationary Euler flow. It is not clear either what is the optimal outflux velocity 
V. From (|A.5P and the discussion just above, it appears desirable to minimize \V\, since 
this minimizes both drag and the imbalance between flow Voo just outside the boundary 
layer and the ambient flow at infinity. On the other hand, we expect that stability becomes 
more delicate in the characteristic limit F — > 0~ , in the sense that the size of the basin of 
attraction of the boundary layer shrinks to zero (recall, we have ignored throughout our 
analysis the size of the basin of attraction, taking perturbations as small as needed without 
keeping track of constants). These would be quite interesting issues for further investigation. 
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